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We theoretically study the (1, 1) triplet to (0,2) singlet relaxation rate in a lateral gate-defined 
double quantum dot tuned to the regime of Pauli spin blockade. We present a detailed derivation 
of the effective phonon density of states for this specific charge transition, keeping track of the con- 
tribution from piezoelectric as well as deformation potential electron-phonon coupling. We further 
investigate two different spin-mixing mechanisms which can couple the triplet and singlet states: a 
magnetic field gradient over the double dot (relevant at low external magnetic field) and spin-orbit 
interaction (relevant at high field), and we also indicate how the two processes could interfere at 
intermediate magnetic field. Finally, we show how to combine all results and evaluate the relaxation 
rate for realistic system parameters. 
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I. INTRODUCTION 

The last decade has seen a great interest in spin 
qubits hosted in semiconductor quantum dots, moti- 
vated by the prospects of easy scalability, weak cou- 
pling to external perturbations, and flexible tunabilityP^ 
Experimental advance has been substantial in the past 
years, and essential operations including qubit initializa- 
tion, manipulation, and readout have been convincingly 
demonstratedP"^This progress is not only exciting in the 
context of quantum computation and information, but it 
also provides a unique platform for studying fundamental 
quantum properties of nanoscale systems. Ongoing effort 
is therefore directed at improving the quality of the spin 
qubits, mainly by trying to reduce qubit dephasingP and 
increase the measurement fidelity?! 

A common method to read out a quantum dot spin 
qubit relies on the so-called Pauli spin blockadeP A dou- 
ble quantum dot is tuned to a (1, 1) charge state, meaning 
that each dot contains exactly one excess electron,^ and 
two of the four resulting (1,1) spin states are used as a 
qubit basis. After qubit manipulation, the double dot po- 
tential is tilted such that a (0, 2) charge state becomes the 
two-electron ground state. For not too strong tilting, the 
only accessible (0, 2) state is a spin singlet, which makes 
the (1,1) — > (0,2) charge relaxation spin-selective. If 
the two qubit basis states contain a different spin-singlet 
component, then one can use charge detection to mea- 
sure the qubit 's final state before tilting, either by do- 
ing transport measurements couplin g th e doubly occu- 
pied dot strongly to an outgoing leadp^ or by detecting 
the charge state with a nearby charge sensorP^H^] 

The accuracy of such a readout depends crucially on 
the effectiveness of the spin blockade: Any leakage out 
of the blocked triplet states reduces the readout visibil- 
ity and thereby distorts the measurement!^ A detailed 
understanding of the spin-flip relaxation responsible for 
such triplet leakage is thus essential in the context of spin 



qubit measurement. Most existing theoretical work along 
these lines was done for single-dot spin relaxations ^ or 
consists of numerical studies of the relaxation rates liZliU 
A thorough analytical study of interdot spin-flip charge 
relaxation in double quantum dots is still missing. 

Here, we study in detail the (1, 1) triplet to (0, 2) sin- 
glet decay rate for a lateral double quantum dot. We 
investigate two spin-mixing mechanisms which provide a 
coupling between the otherwise orthogonal states: (i) For 
small externally applied magnetic fields the coupling is 
believed to be dominated by the effective magnetic field 
gradient over the two dots caused by the hyperfine cou- 
pling of the electron spins to th e ra ndomly fluctuating 
nuclear spins in the host materialpm' (ij) At larger fields, 
this coupling is suppressed for the two polarized triplet 
states for which spin-orbit interaction takes over as dom- 
inating spin-mixing mechanismP^l The second ingredient 
necessary for a finite leakage rate is the dissipation of the 
energy difference A between the initial (1,1) triplet and 
final (0, 2) singlet state, which we assume to be provided 
by the coupling to acoustical phonons in the host mate- 
rial. We derive the function Pi (A) for this specific charge 
transition, which gives the probability that the transition 
is accompanied by the dissipation of energy A by a single 
phonon (alternatively one could call this function the ef- 
fective phonon density of states for the charge transition). 
We include the contribution from piezoelectric as well as 
deformation potential electron-phonon coupling, and we 
find that the piezoelectric contribution to Pi (A) is linear 
at low energies and oc A -5 at high energies, whereas the 
coupling to the deformation potential leads to a contribu- 
tion oc A 3 for low energies and oc A -1 for high energies. 
We finally evaluate explicit relaxation rates using param- 
eters of the experiment of Ref. [T9l We find a decay rate 
~ MHz, which agrees with experimental observations.^ 

Our results are not only relevant for spin qubit readout. 
In Ref. [19] it was suggested that the spin-orbit coupling of 
the (1, 1) triplet and (0, 2) singlet states could also be uti- 
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lized to drive off-resonant microwave-stimulated Raman 
transitions within the (1,1) space. In that case transi- 
tions between the (1, 1) and (0, 2) states would contribute 
to qubit dephasing, and a detailed understanding of the 
mechanisms responsible for these transitions would be 
essential in this context as well. 

The rest of this paper is structured as follows. In Sec.|TT] 
we introduce our model of the double quantum dot and 
present an effective Hamiltonian defining the basis we 
will work in. In Sec. |III| we then investigate the two spin- 
mixing mechanisms (spin-orbit interaction and a mag- 
netic field gradient over the double dot) and we derive 
the matrix elements needed to calculate the leakage rate. 
In Sec. |IV] we study the coupling to the phonon bath 
in detail. We start from the standard Hamiltonian de- 
scribing the electron-phonon coupling and derive from it 
Pi (A) for the (1,1) to (0,2) spin-flip charge transition. 
Finally, in Sec. |Vj we evaluate the leakage rate explicitly 
with realistic experimental parameters. 



II. MODEL 

We consider a lateral double quantum dot, gate-defined 
in a two-dimensional electron gas (2DEG) formed at the 
interface of an AlGaAs/GaAs heterostructure which has 
been grown in the crystallographic (001) direction. The 
setup we have in mind is sketched in Fig. [T^,: The two 
dots (gray circular areas) are separated by a distance d. 
Each dot is approximated to be a two-dimensional har- 




FIG. 1: (a) A lateral double quantum dot tunnel coupled to a 
left and right lead. Two nearby gate electrodes with applied 
voltages Vl,r can change the potential offset of the dots. An 
in-plane external magnetic field B cxt is applied in a direction 
parallel to the double-dot axis, (b) Charge stability diagram 
of a few-electron double dot. We investigate the gray region, 
where both (1, 1) and (0,2) can be stable charge states, (c) 
The spectrum of the electronic Hamiltonian (J2J as a function 
of the detuning e, where we chose t = 0.2 (all energies in units 
of -Bext). The detuning axis e is also indicated in (b). We will 
focus on the regime indicated with the red rectangle, where 
the two lowest-energy states are IS02} and \T + ). 



monic potential well with a typical width a, correspond- 
ing to a single-dot level spacing Hlu — h 2 /2ma 2 where m 
is the effective electron mass (m ~ 0.067m c in GaAs). 
Nearby gate electrodes with applied voltages Vl and Vr 
couple capacitively to the dots and can change the po- 
tential offset of the two dots. We choose the z-axis to be 
parallel to the interdot axis and assume an external mag- 
netic field -B cx t to be applied in the same z-direction (as 
was the case in the experiment of Ref. IT§|) . Two nearby 
leads, labeled L and R, are tunnel coupled to the two 
dots and act thereby as electron reservoirs. 

Part of the few-electron charge stability diagram of 
such a double dot is shown in Fig. [TJd. In the plane 
(Vl,Vr) regions of stable charge configuration {uljUr) 
have a hexagonal shape, where tilir) denotes the number 
of excess electrons on the left (right) dot.^ The regime 
we want to investigate in this work is the shaded trian- 
gle close to the (1, l)-(0, 2) boundary. Here the ground 
state of the system is a (0, 2) charge state, but excited 
(1,1) states are also stable in the sense that they can- 
not decay through sequential tunneling processes like 
(1, 1) — > (0, 1) — > (0, 2). We also assume the tunnel barri- 
ers to the leads to be high enough that the corresponding 
cotunncling processes do not contribute significantly to 
charge relaxation. In that case, we can regard the double 
dot to be a closed system containing two electrons which 
have either a (1, 1) or (0,2) charge configuration. 

The single-dot orbital level spacing hid ~ meV ~ 10 K 
is typically comparable to the charging energy of the 
dot, i.e., the Coulomb energy it costs to add an ex- 
tra electron to the dot. This is such a large energy 
scale that we can disregard all electronic states involv- 
ing higher orbital states and focus on electrons in the 
orbital ground state. Explicitly, the two-dimensional 
ground state wave functions in the left and right potential 
wells read ^(r) = (l/y/2na) exp{-[x 2 + (z + cZ) 2 ]/4a 2 } 
and ^b{v) = (l/^/2na) exp{-[x 2 + z 2 ]/4a 2 }. Using a 
Hund-Mulliken approach we orthonormalize these two 
wave functions, resulting in the basis states^ 

*L(fl)(r) - j , (1) 

with the factor g = (1 — \f\ — s 2 )/ s and the overlap in- 
tegral s = J drtp* L (r)tp R (r) = exp{-d 2 /8a 2 }. 

Including spin into the picture, we can construct from 
these basis states a (1,1) spin-singlet state \S), three 
(1, 1) spin-triplet states \T ) and \T±), and a (0,2) spin- 
singlet state IS02). We assume that in the basis spanned 
by these five states the Hamiltonian describing the ki- 
netic and potential energy of the two electrons as well as 
their Coulomb interaction can be written as^l 

H = H + H t , (2) 

Ho = Bext{|T_) (T_| - \T+) (T + |} - e \S m ) (S 02 \ , (3) 

H t = t{\S){S 02 \ + \S <a ){S\}, (4) 

where we included a Zeeman term coupling to the spin of 
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the two electrons. The field B ext is written in units of en- 
ergy (the sign chosen reflects the fact that the g-factor in 
GaAs is negative, so that a positive B ext corresponds to 
a positive field along the z-axis) and e describes the de- 
tuning between 16*02) and the unpolarized (1, 1) states.^ 
Increasing e can be effected by changing Vl and Vr as 
indicated with the arrow in Fig. [T]d. 

In Fig. [lp we plot the spectrum of H as a function of 
the detuning e, where we chose t = 0.2 (all energies in 
units of B cxt ). Our regime of interest (the gray area in 
Fig.]!^) is where 15*02) is the ground state, indicated with 
the red dotted rectangle in Fig. [TJi. We assume that in 
this regime e 3> t, so that the tunneling Hamiltonian H t 
can be treated as a perturbation. An excited (1, 1) state 
can only decay to the (0, 2) ground state if its wave func- 
tion is a spin-singlet (or contains a singlet component). If 
the system is in a pure (1,1) spin-triplet, it cannot decay 
and stays blocked in the excited state. 

We now focus on such a spin-blocked situation and as- 
sume that the system is initially in \T + ). The purpose 
is to calculate the relaxation rate from \T + ) to 15*02) (see 
the wiggly line in Fig. [TJ:). For this we need two in- 
gredients, which we will investigate in detail in the next 
two Sections: (i) We need a perturbation H sm with a fi- 
nite matrix element between the states \T + ) and |5o 2 ). 
(ii) The energy difference between initial and final state 
Et + — Es 02 = A has to be dissipated by the environment 
of the double dot, for which we assume the coupling to 
acoustic phonons to be responsible. 



III. SPIN-MIXING PERTURBATIONS 

As noted before, \T + ) and IS02) are orthogonal in spin 
space, and the perturbation _ff sm thus has to be of a 
spin-mixing nature. We will investigate two such per- 
turbations: (i) spin-orbit interaction, which mixes the 
spin and orbital parts of the electrons' wave functions, 
and (ii) a magnetic field gradient over the double dot, 
i.e., a difference between the effective magnetic fields at 
the positions of the left and right dot. Below we will in- 
troduce the two perturbations and calculate the resulting 
matrix elements coupling |T+) and |So 2 )- 

A. Spin-orbit interaction 

Spin-orbit interaction perturbs single-particle states in 
the two dots resulting in mixed spin-orbital eigenstates 
instead of pure spin statesP Thereby it can give rise to 
"spin-fli p" tunne l coupling of states with apparent oppo- 
site spinP^M25j p or eacn electron the spin-orbit Hamil- 
tonian reads 

H so = a(-pyfr x +Pz&y) +P(-p3:Vx +Py&y), (5) 

where p is the momentum of the electron and & x y z are 
the three Pauli matrices. The first term in |5]) is the 
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FIG. 2: The orientation of our zz-plane with respect to the 
(100) and (010) crystallographic directions. 

so-called Rashba term and the second the Dresselhaus 
term. For typical 2DEG's in GaAs the corresponding 
spin-orbit length l so , i.e. the distance an electron has to 
travel to have its spin rotated by ~ 1, is of the order 
l so ~10 /xm, usually much larger than the size of the 
dotsP^H The ratio of the two parameters, a//3, depends 
on the detailed confining potential of the 2DEG and can 
in practice be smaller as well as larger than 1. 

The spin-orbit Hamiltonian (|5]) is written such that 
the x-, y-, and z-axes point respectively along the (100), 
(010), and (001) crystallographic axes, and it is assumed 
that the 2DEG lies in the ij/-plane. Transforming this 
Hamiltonian to the coordinate system of Fig. [I] we find 

H so = a[-p x a z + p z a x ] 

+ j3 [p z (a x sin 2% - a z cos 2%) 

+ p x (a x cos2x + a z sm2x)}, (6) 

where x is the angle between the double-dot axis and the 
(100) crystallographic direction, see Fig. [2} 

In the two-electron position representation, the basis 
states |T + ) and |5o 2 ) rea d explicitly 

(ri )2 |T+) = ^{*it(ri)*jtr(r 2 ) - ^ t (ri)* it (r 2 )}, 
(1X2I602) = ^{*Rt(ri)*iy.(r 2 ) - *i?;(ri)-Ffl T (r 2 )}, 

where the vE'(r) are now in fact two-component spinors. 
Using that p = —ihd r we find straightforwardly that 
only terms with p z have a non- vanishing matrix element 
between |T + ) and IS02), the total matrix element reading 

T so = (S 02 \H ao \T + ) = -ihJL J^-s ia + psm2 X ). 

(7) 

We see that T so depends non-trivially on the angle x- For 
a double dot with its interdot axis pointing in the crystal 
(110) direction the two spin-orbit terms add construc- 
tively and T so oc (a + whereas dots with the interdot 
axis along the (110) direction have T so oc (a — f$). 

We also note here that the direction of the external 
field B ext can have a great influence on the effectiveness 
of the coupling. For instance, for a double dot with the 
interdot axis constructed along the (110) direction (i.e., 
X = tt/4), the Hamiltonian ([6]) reads 

H so = a(-p x a z +p z a x ) + /3(p z a x +p x a z ). (8) 
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If the external field now points in the x-direction, then 
all terms with p z come with a x , which cannot provide a 
spin-flip. In this case we have thus T so = 0. 

There is however another, second-order, spin-orbit me- 
diated process which can take place: H so couples the 
ground state \T + ) to a (1,1) spin-singlet which involves 
an excited orbital state. This excited state can be as- 
sumed to be coupled to IS'02) with a coupling energy of 
~ t. The spin-flip now takes place inside one of the dots, 
and does not depend on the orientation of the interdot 
axis. An estimate of the magnitude of the resulting ef- 
fective matrix element coupling \T + ) to 16*02) gives 

T so .2~ (a + P) k -^, (9) 
a nui 

If we assume that t ~ p,eV and tko ~ meV, then a com- 
parison with the T so found above yields that the "direct" 
coupling dominates as long as d/a < 8 is the case.^ Of 
course, if T so = (due to a special orientation of the in- 
terdot axis and B ext ) then T SOj 2 still provides a spin-orbit 
assisted spin-mixing coupling between \T + ) and IS'02) • 

B. Magnetic field gradient 

A magnetic field gradient over the two dots generally 
mixes all four (1,1) states. All triplet states acquire a 
spin singlet component, and thus are coupled to IS'02) by 
the tunneling Hamiltonian H t . Such gradients coul d be 
due to a deliberately fabricated on-chip micromagnetP^D 
or to the effective magnetic fields set up by th e ran domly 
fluctuating nuclear spins of the host material) 10 * 28 * 

The Hamiltonian describing the coupling of two dif- 
ferent effective magnetic fields SBl and <5B^ to the two 
electrons reads 

H gI = - 5B L ■ S L - SB R ■ S R , (10) 

where we again have set gp R = — 1 for convenience, and 
Sl(r) is the (dimensionless) spin operator for the electron 
in the left(right) dot. In the basis of the (1, 1) singlet and 
triplet states this Hamiltonian readiP^ 

^ = E {-^f i T °> < T ±i ± 7! 1 5 ) < r ±i + H - c j 

-6B',{\T + ){T+\-\T-){T-\} 

-6B*{\S)(T \ + \T ){S\}, (11) 

where we used the notation SB^ a = SB^ a ± iSB^ a . The 
symmetric and antisymmetric fields we used are defined 
as SB S = \{SB L + SB R ) and SB a = \{5B L - 5B R ). 

We assume that the fields 5B^ and SB R are much 
smaller than the externally applied field. Then we can 
use first-order perturbation theory to find the singlet 
admixture in |T + ) caused by H gr , which yields |T + ) ss 
\T+) - {SB+/V2B cxt )\S). Therefore, the spin-flip ma- 
trix clement due to the perturbation H t + H gL - reads 

8B+ 

T SI = (S 02 \H t + H er \T + ) = -t-^-. (12) 



If the gradients are caused by randomly fluctuating nu- 
clear spins, the typical magnitude of the effective fields is 
approximately SB^ r ~ A/\N, where A is the material- 
specific hypcrfinc coupling energy and is the number 
of nuclei in each dot. For GaAs dots A ~ 100 /ieV and 
typically N ~ 10 6 , implying that 5Bl r is in the regime 
of 1-5 mT, which has been confirmed experimentally.^ 



IV. ELECTRON-PHONON COUPLING: 
SINGLE-PHONON PROBABILITY 

Both perturbations outlined above provide effectively 
a coupling between \T + ) and IS'02) and can thus cause 
relaxation. The only ingredient still missing to fully de- 
scribe transitions between the two levels, is a mechanism 
dissipating the energy difference A between initial and fi- 
nal state, typically 50-500 /ieV. We assume that this en- 
ergy is absorbed by the (acoustical) phonon bath in the 
host material. Often the contribution to the electron- 
phonon coupling from the deformation potential is ne- 
glected, which is generally justified for phonons with en- 
ergies below ~ 10 KP3 However, phonons with a wave 
vector larger than the inverse in-plane dot size 1/a are 
emitted almost exclusively in the y-direction (perpendic- 
ular to the plane of the 2DEG), and emission of piezoelec- 
tric phonons in this dire ction is strongly suppressed by 
the crystalline anisotropy! 14 * 15 * Since for a typical GaAs- 
hosted double quantum dot system a phonon wave vector 
of 1/a corresponds to an energy of ~ 100 fJ.eV, we will 
keep in our calculation both the coupling to piezoelectric 
phonons and deformation phonons. 

The Hamiltonian describing the coupling between the 
electrons (density operator p) and phonons (creation and 
annihilation operators fi^) reads 

#c- P h = ^q,p^q[^q,p + °-qjJ' 

where A qiP are the coupling matrix elements and /5 q — 
J dr e~ 4q r /5(r) is the Fourier transform of the electronic 
density operator. The sum runs over all allowed phonon 
wave vectors q and includes three polarizations (one lon- 
gitudinal and two transversal), labeled by p = l,t%,t2- 
We will neglect the mismatch of phonon velocities at the 
GaAs-AlGaAs interface, and treat the phonon bath as 
that of bulk GaAs. 

The matrix elements for electron-phonon coupling read 

where p is the mass density (p = 5.3 x 10 3 kg/m 3 for 
GaAs), V the normalization volume, and we will assume 
that the phonons have an isotropic linear dispersion re- 
lation at all energies of interest, i.e., w qiP = v p q, with v p 
the polarization-dependent sound velocity. 
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The constant M 



(p) 



M, 



(p) 

pe 



M(j c f contains a contri- 



bution from both types of electron-phonon coupling, 
AfW = ie/ii 4 A qiP , 
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M def = Sg(5 p ,;, 



(15) 
(16) 



where the coupling to the deformation potential only in- 
volves longitudinal phonons, expressed by the (5-function 
in (16). We used here the piezo-electric constant, hi4 — 



1.38 x 10 9 V/m in GaAs, and the deformation poten- 
tial, which is 5 = 13.7 eV for GaAsP 3 The coupling to 
piezoelectric phonons involves the anisotropy factors 



A — 

q,p ? 



qzqy 



(p) , (p) 
qzqxe'y + qyQze- 



(17) 



where e( p ) is the unit polarization vector for the polariza- 
tion p. The factors as written in ( 17 ) are in the coordinate 



system of the crystal structure, i.e. the i-direction along 
(100), y along (010), and ~z along (001). We would like 
to relate the factors to the coordinate system of Fig. [TJi. 
We thus write in terms of the spherical coordinates of q 



Aqj = 9 cos 2 (6») sin 4 (6») sin 2 (20 + 2 X ) 
A q>t2 = sin 2 (2(9) cos 2 (20 + 2 X ) 



(18) 



qtl = i[l + 3cos(26»)] 2 sin 2 (6»)sin 2 (20 + 2x), (19) 



(20) 



where 9 = corresponds to q parallel to our y-axis, and 
<f> gives the azimuthal angle of q in our zz-plane. The 
angle x is the angle between the double dot axis and the 
crystallographic (100) direction: A wave vector q with 
given cj> thus has an azimuthal angle (f> + X m the crystal's 
coordinate system (see Fig. [2]). 
We see that we can write 

|A q , P | 2 = (9$A«,p + W-W) , (21) 

with the two dimensionless coupling constants, 
_ (e/> 14 ) 2 



g(p) 



and g de f = 



2ir 2 hpvf a 2 ' 



(22) 



the latter being dependent on the dot size a. 

The relaxation rate T of the excited (1,1) triplet state 
to the (0, 2) ground state will be calculated using a second 
order Fermi's golden rule, 



2tt 

T 



V (f\H'\v)(v\H'\i) 



5{E f - Ei), (23) 



where H' = Hc^+Hsm, with H sm being one (or both) of 
the spin-mixing Hamiltonians presented in Sec. |III| The 
energy difference A between initial state \T + ) and hnal 
state ISba) (which is equal to the energy of the emitted 
phonon) is assumed much larger than the temperature 
and we therefore take as initial state a direct product 
of |T_|_) and the phonon vacuum |vac), and as final state 
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|lqp)i where one phonon with wave vector q and 



polarization p has been created. 

From the explicit wave functions of |T+) and 15*02 
we calculate the diagonal matrix elements (T + \H, 
and (S 2\H e . ph \S 02 ), and fincP 1 



e-ph|2 + ) 



q./j 



2Fq 

A 



iq-d^* 



q,p 



-A), 
(24) 

where the spin-mixing matrix element T = {So2\H sm \T + ) 
and we used the Fourier transform of the squared elec- 
tronic wave function F q = J dr e~ 4q ' r |<I'ft(r)| 2 . We can 
evaluate this Fourier transform explicitly and, anticipat- 
ing that the (5-function enforces hu)q,p = A, we write 

r - 2?r 4 V |T|2 i,\ \ 2 e~- 2 \llW.) 

q,p 

x sin (^q z d) 2 S(hujq t p — A). (25) 

The exponential function exp{— a 2 (q 2 + q 2 )} suppresses 
the contribution from phonons having a wave vector with 
in-plane components larger than the inverse system size 
1/a. Indeed, the electronic density profile P q is expo- 
nentially small for these wave vectors. The sine function 
s'm(^q z d) 2 describes the interference between the cou- 
pling to an electron in the left and right dot: A phonon 
wave with given wave vector q has a phase difference 
S4> = qzd between the two dot positions.^ 

We convert the sum over q into an integral and then 
finally find that we can write for the relaxation rate 



2tt 



T| 2 Pi(A), 



(26) 



where the function Pi (A) gives the total probability that 
the energy A is absorbed by a single phonon, either by 
piezoelectric coupling or by coupling to the deformation 
potential. Alternatively, one could call this function the 
effective phonon density of states for the (1,1) triplet to 
(0, 2) singlet charge transition. 

The total single-phonon probability is the sum of the 
contributions from the different types of phonons, 



P 1 (A) = P Mcf (A) + ^P 1 (p p ) e (A). 



(27) 



For the piezoelectric phonons we find 



irSe(A) = 



i P J r /2 ,„ sin6» „ /Adsinfl 



gpe 

A 



d6 



2„2 



x exp 



A 2 a 

h 2 v 2 



sin 2 6 } , (28) 
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where the dimensionless functions f p (x) are 
ft(x) = ^cos 2 (0)sin 4 (%_(x), 



(29) 



f n (x) = -[l + 3cos(2f?)] 2 sin 2 (%_(x), (30) 



/ t2 (x) = -sin 2 (2%+(x), 



(31) 



in terms of the function 



// (.<■) = j -./„(./■) ( ? 2 - ; cos(4\ )./„(.!■ i 



± 



cos(4x)Ji(as), 



(32) 



with J„ (x) the n-th order Bessel function of the first kind. 
The contribution from the coupling to the deformation 
potential reads similarly 



Pl,def(A) 



5dcf 

IT 



i 



exp 



/ Ad sin 9 




The remaining integral over the polar angle 9 has to 
be evaluated numerically. We can however arrive at an- 
alytical results in the limits of small and large phonon 
energies. For small energies, meaning Tiv/A ^> a,d, the 
single-phonon probabilities Pi can be expanded in A. 
Setting v t i ~ Vt2 = Vt we then find to leading order 



p£U A ) = 

Pl,def(A) = 



G 



(0 
3po 



105(1 - s 2 ) A 

° 9pe 



dA 



dAY 



105(1 - s 2 ) A \%v t J ' 
In; 



g dci ( a A \ 2 / dA N ~ 



6(1 -s 2 ) A 



(34) 
(35) 
(36) 



where P 



(t) 



P 



(ti) 



P 



(t2) 



i, pc — 1 l.pc t J i.pc- We find for small energies 
a linear piezoelectric Pi. po (A) and a cubic deformation 
Pi t dcf{A), meaning that the phonon bath is superohmic 
in this setup. The result for the piezoelectric phonons 
agrees up to a prefactor with previous calculations of the 
phonon density of states for the (1,0) to (0,1) charge 
transition where all anisotropy factors were set to onepa 
In the opposite limit of large energies, hv/A <C a, we 
find qualitatively different results compared to Ref. 
We see from the exponential factors in (p8| and (33) 



that in this regime only very small angles 9 are rele- 
vant. Indeed, in this case only the confinement in the 
^-direction is strong enough to create an electronic den- 
sity profile with non- vanishing Fourier components of the 
order ~ A/hv, and phonon emission takes place almost 
exclusively in the y-direction. For the piezoelectric cou- 
pling the anisotropy factors ^4q lP now bring in small fac- 
tors of which cannot be ignored (see also Ref. ITlj) . Since 



in this limit Ai oc 9 4 and A t i,A t 2 oc 6 2 , we expect the 
dominating piezoelectric contribution for large energies 
to come from transversal phonons. 

To evaluate the single-boson probabilities in this limit, 
we expand sin 9 sa 9 and extend the range of integration 
over 9 from to oo. Then we find that to leading order 



PSe(A) 



Pl,def(A) 



2a 2 (1 - s 2 ) 



1$ 




A 





ffdcf 

2A ' 



(37) 
(38) 



The contribution from longitudinal piezoelectric phonons 
is Pij ((7po / A)(hvi/aA) e , which is much smaller than 
Pi pe an d therefore ignored. The large-energy result ( 37 ) 



for the piezoelectric phonons is qualitatively different 
from the results presented in Ref. [29l which predicted 
that Pi (A) (g pe / A)(Tiv /aA) 2 , the difference arising 
from the inclusion of the anisotropy factors. 



We see from ( |37[ ) and ( 38 ) that the contribution from 
deformation phonons becomes important when A^/a ~ 




FIG. 3: The total single-phonon probability Pi = Pi, pc + 
Pi.dcf as a function of A. The phonon energy A is plotted in 
units of A a = hvt/a and Pi in units of P a = a(eh\4,) 2 /~h 2 vf p. 
For all plots we have set vi/v t — 1.73 and E/ae/114 = 0.50. 
(a) The total Pi (A) at x — f° r three different size ratios 
d/a. (inset) The dependence of the maximum Pi for d/a = 5 
on the angle \. (b) The total Pi (A), as well as the separate 
contributions from piezoelectric and deformation phonons, for 
X = and d/a = 5 on logarithmic scales. The expected power 
laws are included as guides to the eye. 
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1/2 



hvt ■^vfehi_4 l /v^'S, which is approximately 1.4 meV-nm 
using realistic parameters for GaAs!^ For a dot size of 
a = 20 nm (an orbital level spacing of « 1.4 meV) 
we find that the relevant energy scale is A ~ 300 peV, 
which indeed lies inside our regime of interest. 

In Fig. 3] we plot the single-phonon probability Pi (A) 
for typical double-dot parameters. In all plots A is renor- 
malized to units of A = hvt / a and Pi is plotted in units 
of P a = a(eh 1A ) 2 /h 2 v?p. For v t = 3.0 x 10 3 m/s and 
a = 20 nm we find A a = 99 peV. The parameter vi/v t 
was set to 1.73 and the ratio E/aehu to 0.50. In Fig. [3^, 
we show the total Pi at x = for three different size ra- 
tios d/a. The maximum of Pi always occurs on the scale 
A ~ A a , where the wave length of the emitted phonon is 
comparable to the system's in-plane dimensions. At low 
energies Pi is approximately linear and at high energies 
it is suppressed, ultimately being dominated by the de- 
formation contribution making Pi cx A -1 . In Fig. we 
plot the total Pi as well as the two separate contribu- 
tions on logarithmic scales (for x = and d/a = 5), and 
we added guides to the eye corresponding to the power 
laws expected in the different limits. The blue dotted 
line shows the piezoelectric contribution. Up to a few 
A a this contribution indeed dominates, being linear at 
very small energies. For A > A a it becomes suppressed 
as oc A~ 5 and at higher energies the dominating contri- 
bution comes from the coupling to the deformation po- 
tential, the green dotted line (see also Eqs. 37 and 38 1. 
The inset of Fig. [3^, shows the dependence of the max- 
imum of Pi on the angle x (for d/a = 5). We see that 
the density of states depends on the relative orientation 
of the double dot axis and the (100) crystal direction, its 
variation being however only ~ 1 %. 



(a) 



.1 4 




50 100 150 200 




A = 50 ,u,eV 



FIG. 4: Relaxation time Trei from |T+) to IS02) when relax- 
ation is mediated by spin-orbit interaction. To make this 
plot we used vi = 5.2 x 10 3 m/s, v t = 3.0 x 10 3 m/s, 
hu = 1.38 x 10 9 V/m, p = 5.3 x 10 3 kg/m 3 , x = t/4, and 
a = P = 100 m/s. (a) T Tcl for a = 20 nm and d = 100 nm. 
(inset) The relaxation rate V for the same parameters calcu- 
lated at A = 50 peV as a function of \. (b) T rc i as a function 
of d/a at two different energies for fixed a — 20 nm (solid 
lines) and d — 150 nm (dashed lines). 



V. SPIN-FLIP CHARGE RELAXATION RATE 



Now we can use Eqs. (26)-(32) to evaluate explicit re- 
laxation rates for the T + — > S02 transition. For a specific 
experimental setup one can estimate the relative magni- 
tude of the matrix elements given in Eqs. Q and Q, 
and decide which process dominates. Here, we will focus 
on the case of a large external magnetic field, such as 
was the case in the experiments of Ref. [THl We assume 
that P ext is large enough so that |T so | 3> |T gr |. In that 
case the relaxation from \T + ) to IS'02) is mainly caused 
by spin-orbit interaction. 

We thus use T so in (26) and take typical material 
parameters for GaAs. For the dot size we take again 
a = 20 nm and the interdot distance is set five times as 
large, d = 100 nm. In the experiment of Ref.H~9lthe inter- 
dot axis was fabricated along the crystal (110) direction, 
so we set x — tt/4. For this angle the two spin-orbit terms 
add constructively, and we choose a = /3 = 100 m/s, 
such that l so = h/m(a + /3) « 8.6 pm. In Fig. [4^, we plot 
the resulting relaxation time T rc i = T -1 , which is found 
to be typically T rc i ~ 1 /is, the order of magnitude of 
which agrees with experimental observations PB The in- 



set to Fig. [4^, shows the dependence of the relaxation rate 
r on the angle x at A = 50 ^eV, close to the minimum 
relaxation time. We see that the rate indeed vanishes 
for the angles x — 37r/4, 77r/4, where the Rashba and 
Dresselhaus terms add destructively. 

In Fig. [4]d we show how the relaxation time depends on 
the size ratio d/a, for two different energies A = 50 /xeV 
(red lines) and A = 250 ^eV (blue lines). The dashed 
lines have a fixed interdot distance d = 150 nm and 
the solid lines have a fixed dot radius a — 20 nm. All 
other parameters are the same as in Fig. [4^,. We see 
that a large size ratio d/a suppresses the relaxation effi- 
ciently: For widely separated dots the overlap of the two 
single-dot wave functions becomes exponentially small, 
and this suppresses the matrix element T so . In the limit 
of strongly overlapping wave functions, i.e. d/a going to- 
wards 1, we see that relaxation is much more efficient 
for the smaller system with d = a = 20 nm. Indeed, for 
d = a = 150 nm we find A a s» 13 peV, so in this case both 
energies are larger than A where the electron-phonon 
coupling matrix elements are suppressed. 



For smaller external magnetic fields, or other angles \i 
one could be in the situation where the spin-orbit and 
field-gradient give rise to matrix elements of the same 
order of magnitude, |T so | ~ |7g r |. In this case one has 



to use in ( 26 ) the total matrix element T ti 



gi'- 



which possibly includes interference terms between the 
two mechanisms, 



71, 



-iSBl + SBl hds(a + (3 sin2 X ) \ 
V2B cxt 4a 2 VT^ ' ' 



We see that the spin-orbit mechanism interferes with the 
y-component of the difference field <5B a . By tuning SB^ 
or B ext one could thus enhance or counteract the spin-flip 
tunneling enabled by spin-orbit interaction. One word of 
caution is however required here: If the field gradients 
are caused by the effective hyperfine fields, then the fi- 
nal state |/) in (23) is different for a spin-orbit and a 
hyperfine mediated transition. Indeed, in the course of 
hyperfine induced spin-flip tunneling the spin of one of 
the nuclei is raised by h, which does not happen during 
a spin-orbit mediated transition. In that case one has to 
calculate separately the two contributions to the relax- 
ation rate @ or, equivalently, use \T\ 2 = \T so \ 2 + |T gr | 2 . 



dot tuned to the Pauli spin blockade regime. We first de- 
rived an effective phonon density of states -Pi (A) for this 
charge transition, and found that at small energies Pi is 
linear in energy, cx A, and dominated by the piezoelec- 
tric electron-phonon coupling, whereas at large energies 
the Pi is dominated by the coupling to the deformation 
potential and is cx A -1 . Then, we investigated two differ- 
ent spin-mixing mechanisms coupling the spin triplet and 
singlet states: a magnetic field gradient over the double 
dot (relevant at low external magnetic field) and spin- 
orbit interaction (relevant at high field) . We showed how 
the spin-orbit-mediated coupling depends on the device 
geometry as well as on the in-plane direction of the ap- 
plied magnetic field. We finally combined all results and 
took realistic system parameters to evaluate the explicit 
detuning-dependent relaxation rate, which we found to 
be of the order of ~ MHz. 
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